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Abstract: The relation between isoenergetic and Hamil- 
tonian thermostats is studied and their equivalence in 
the thermodynamic limit is proved in space dimension 



d = 


1,2. 
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I. THERMOSTATS 

In a recent paper [1] equivalence between isokinetic 
and Hamiltonian thermostats has been discussed hcuris- 
tically, leaving aside several difficulties on the under- 
standing of the classical dynamics of systems of infinitely 
many particles. Understanding it is, however, a necessary 
prerequisite, because strict equivalence can be expected 
to hold only in the thermodynamic limit. In this paper 
we proceed along the same lines, comparing the isoen- 
ergetic and the Hamiltonian thermostats, and study the 
conjectures corresponding to the ones formulated in [1] 
for isokinetic thermostats, obtaining a complete proof of 
equivalence in 1 and 2-dimensional systems. 

Here the class of models to which our main result ap- 
plies is described in detail. The main result is informally 
quoted at the end of Sec. I after discussing the physics and 
the equations of motion of the models; a precise state- 
ment will be theorem 1 in Sec. Ill and it will rely on a 
property that we shall call local dynamics: the proof is 
achieved by showing that in the models considered the lo- 
cal dynamics property holds as a consequence of the the- 
orems 2-9, each of which is interesting on its own right, 
discussed in the sections following Sec. III. 

A classical model for noncquilibrium statistical me- 
chanics, e.g. see [2], is a test system in a container Oo, 
and one or more containers Clj adjacent to it and enclos- 
ing the interaction systems. 



cV 2 



c 3 



x — (Xq, Xq, Xi, Xi, . . . , X^, X„) 



Fig.l: The 1 + v boxes Q d n A, j 
Co, Ci, ■ ■ ■ , C v and contain No , Ni , . . . 
Hons and velocities denoted Xo,Xi, . 
X„, respectively. 



= 0, . . . , v, are marked 
, N„ particles with posi- 
. , X„, and X , Xi, . . . , 



To fix the ideas the geometry that will be considered 
can be imagined (see Fig.l, keeping in mind that it is 
just an example for convenience of exposition and which 
could be widely changed) as follows: 

(1) The test system consists of particles enclosed in a 
sphere ilo = S(Do) of radius Do centered at the origin. 

(2) The interaction systems consist of particles enclosed 
in regions Clj which are disjoint sectors in M. d , i.e. disjoint 
semiinfinite "spherically truncated" cones adjacent to f^o, 
of opening angle u)j and axis k.,: fij ■ = {£ e R d , |£| > 
£>o,£-kj < j = \,...,v ■ 

The initial configurations x of positions and veloci- 
ties will be supposed to contain finitely many particles 
in each unit cube. Thus the test system will consist of 
finitely many particles, while the interaction systems are 
infinitely extended. 

The motion starting from x must be defined by first 
regularizing the equations of motion (which are infinitely 
many and therefore a "solution" has to be shown to ex- 
ist). The regularization considered here will be that only 
the (finitely many) particles of the initial data x inside 
an artificial finite ball A = S(r) of radius r > Do will be 
supposed moving. 

I.e. for the same initial data x only the particles in 
Qq, Q\ n A, . . . , £l v n A will move and kept inside A by an 
elastic reflection boundary condition at the boundary of 
A, while they will never reach the boundaries of the ttj 's 
because of the action of a force, modeling the walls of the 
Qj and diverging near them. 

The particles of x located outside the container Do U 
Uj>o(^j H A) are imagined immobile in the initial posi- 
tions and influence the moving particles only through the 
force that the ones of them close enough to the boundary 
of A exercise on the particles inside A. 

In the "thermodynamic limit 7 , which will be of central 
interest here, the ball A grows to oo and the particles that 
eventually become internal to A start moving: in other 
words we approximate the infinite volume dynamics with 
a finite volume one, called A-regularized, and then take 
an infinite volume limit. 

A configuration x will be imagined to consist of a 
configuration (Xo,Xo) with Xo contained in the sphere 
£(£>o), delimiting the container flo of the test systems, 
and by n configurations (X 
j = 1, 



*■ 3 ' X J ) 



with Xj c Oj n R d ), 



I: Thermostats 



2 



Phase space: Phase space H is the collection of locally 
finite particle configurations x — (. . . ,Qi,qi, . . 

ar=(Xo,Xo,Xi,Xi,... ) X w ,X I ,) = / (X,X) (1.1) 

withXj C fij andqi G R d : in every baWS(r') of radius r' 
and center at the origin O, fall a finite number of points 
of X. //A = OnS(r) and x G TL we shall denote x\ the 
positions and velocities of the particles of x located in A 
together with the positions of the particles of x outside A. 

The particles of x located outside A will be regarded 
as immobile. The particles are supposed to interact with 
each other via a potential p and with the walls with a 
potential ip: 

Interaction: Interparticle interaction will be through a 
pair potential ip with finite range r v smooth, decreasing 
and positive at the origin. The walls of the containers Qj 
are represented by a smoooth decreasing potential ip > 
of range <C r v and diverging as an inverse power of 
the distance to the walls. 

Hence the potential tp is superstable in the sense of 
[3]: a property that will play an important role in the 
following. The value of the potential tp at midrange will 

_ def 

be denoted Tp and < Tp < tp = ip(0); the wall potential 
at a point q at distance r from a wall will be supposed 
to be given by 

1>(q) = (£) a <p , r< r f (1.2) 

with a > and r equal to the distance of q to the wall; 
for larger r it continues, smoothly decreasing, reaching 
the value at r = r^. The choice of ip as proportional 
to ipa limits the number of dimensional parameters, but 
it could be made general. The restriction <C r v is 
required to facilitate the interaction between particles in 
fio and particles in Uj^Ctj. 

The particles in Qq are supposed to interact with all 
the others but the particles in flj interact only with the 
ones in Clj U flo- the test system in f2 interacts with all 
thermostats but each thermostat interacts only with the 
system, see Fig.l. 

The equations of the A-regularized motion (see Fig.l), 
aside from the reflecting boundary condition on the ar- 
tificial boundary of A, concern only the particles in 
Oo U Uj>o(Oj A) and will be 

mX 04 - - diUoi^o) -^9 i [/o,i(Xo,X j ) + * i (X ) 

j>0 

mXji = - diUjiXj) - diUo,jQto,'X.j) - aa 3 ltji (1.3) 

where (1) the parameter a will be a = 1 or a = depend- 
ing on the model considered; 

(2) the potential energies Uj(Xj), j > and, respec- 
tively, f/o,j(Xo,Xj) denote the internal energies of the 



various systems and the potential energy of interac- 
tion between the system and the thermostats; hence for 
Xj C Qj H A the Uj 's are: 

U j (X j )= E 

(3) the first label in Eq.(1.3), j = or j = 1, . . . , v re- 
spectively, refers to the test system or to a thermostat, 
while the second indicates the derivatives with respect to 
the coordinates of the points in the corresponding con- 
tainer. Hence the labels i in the subscripts (j, i) have Nj 
values and each i corresponds (to simplify the notations) 
to d components; 

(4) the multipliers ctj are, for j = 1, . . . , v, 

def Qi 

a i — , tz ; — ^ , x ; — i W1 th 

dN^T^/m' (L5) 

Qj — — Xj • djU j(X , Xj), 

where f Nj k B Tj(xf= K^Xj) d = f f X|and a, are cho- 
sen so that Kj t \ptj) + Uj^Aptj) = Ej t \ are exact con- 
stants of motion if a = 1 : the subscript A will be omitted 
unless really necessary. A more general model to which 
the analysis that follows also applies is in [4]. 

(5) The forces #(Xo) are, positional, nonconservative, 
smooth "stirring forces", possibly absent. 

(6) In the case of A-regularized thcrmostatted dynamics 
we shall consider only initial data x for which the ki- 
netic energies if^A (Xj ) of the particles in the IX, n A's 
are > for all large enough A. Then the time evolu- 
tion is well defined for t < t\(x) where t\(x) is defined 
as the maximum time for which the kinetic energies re- 
main positive (hence the equations of motion remain well 
defined because the denominators in the aj stay posi- 
tive). It will be important to remark that if t\(x) < +oo 
the moving particles positions and velocities have a limit 
even as t — > t\(x) because the accelerations ajqji remain 
bounded even though aj — > oo (by the Schwartz inequal- 
ity a bound on ajqji could be N% | max \d<p\ if Na is the 
number of particles in A). 

The equations of motion with a = 1 will be called 
A-regularized isoenergetically thcrmostatted because the 
energies Ej = Kj + Uj stay exactly constant for j > and 
equal to their initial values Ej . The equations with a = 
in Eq.(1.3) will be considered together with the above and 
called the A-regularized Hamiltonian equations. 

The qualifier "Hamiltonian" refers to case a = in 
which no dissipation occurs even though, strictly speak- 
ing, the equations, unless # = 0, are still not Hamiltonian 
(in spite of ctj = 0). 
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Remark that Qj is the work done, per unit time, by the 
test system on the particles in the j-th thermostat. 

The essential physical requirement that the thermo- 
stats should have a well defined temperature and density 
is satisfied by an appropriate selection of the inital con- 
ditions. The guiding idea is that the thermostats should 
be so large that the energy that the test system trans- 
fers to them, per unit time in the form of work Qj, is 
acquired without changing, not at least in the thermo- 
dynamic limit, the average values of the densities and 
kinetic energies {i.e. temperatures) of the thermostats in 
any finite observation time > 0.. 

To impose the latter requirement, in the thermody- 
namic limit, the values Nj,Ej will be such that [^j^a] 



• ej : with Sj , ej > 



-j^ Sj and, for j > 0, —_ 
fixed in a sense that is specified by a choice of the initial 
data that will be studied, and whose physical meaning is 
that of imposing the values of density and temperature 
in the termostats, for j > 0. 

Initial data: The probability distribution /in for the ran- 
dom choice of initial data will be, if dx = IlJ=o ~~W~T~ 2 ~ > 
the limit as An — > oo of the finite volume grand canonical 
distributions on H 

Ho M {dx) = const e ~ Ho - A °W dx, with (1.6) 
H ,Ao(x) = ^2Pj(K jM (x) - XjN jM + U jM (x)) 

1=0 

fa d M _L_ > o, \j e k, 



Remarks: (a) The values /3n = 



> 0, An € 



also 



fixed, although they bear no particular physical meaning 
because the test system is kept finite. 

(b) Here A = (A , Ai, . . . K) and T = (T , T u . . . , T v ) are 
fixed chemical potentials and temperatures, and An is a 
ball centered at the origin and of radius rn. 

(c) The distribution /in is a Gibbs distribution obtained 
by taking the "thermodynamic limit" An — > oo. Notice 
that is a product of independent Gibbs distributions 
because H does not contain the interaction potentials 

(d) Ao should not be confused with the regularization 
sphere A: it is introduced here and made, right away, oo 
only to define /xn- 

(e) The theory of the thermodynamic limit implies the 
existence of the limit distribution /in, either at low den- 
sity and high temperature or on subsequences, [5]. In 
the second case (occurring when there are phase transi- 
tions at the chosen values of the thermostats parameters) 
boundary conditions have to be imposed that imply that 
the thermostats are in a pure phase: for simplicity such 
exceptional cases will not be considered; this will be re- 
ferred to as a "no-phase transitions" restriction. 

Main result: In the thermodynamic limit, the ther- 



mostatted evolution, within any prefixed time interval 
[0,0], becomes the Hamiltonian evolution at least on a 
set of configurations which have probability 1 with respect 
to the initial distribution /in, in spite of the non station- 
arity of the latter. 
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The initial data will be naturally chosen at random 
with respect to /in. Let the "pressure" in the j-th thermo- 
stat be defined by pj([3, A; A ) = ^-pjnA^I lo § Z i- A <> A ) 
with 



N=0 



dXN 



(2.1) 



-0( — XN+Kj (x N ) + Uj (w)) 



where the integration is over positions and velocities of 
the particles in A n flj. Defining p(f3,X) as the ther- 
modynamic limit, An — > oo, of Pj(f3, A;A ) we shall say 
that the thermostats have densities 5j , temperatures Tj , 
energy densities ej and potential energy densities uj , for 
j > 0, given by equilibrium themodynamics, i.e.: 



d P (l3j,Xj) ! 

v ' »x, ■ >i 



(2.2) 



i = dfJj u j = e i- 2 5 3Pj 

which are the relations linking density Sj, temperature 
Tj = (ksflj) , energy density ej and potential energy 
density Uj in a grand canonical ensemble. 

In general the A-regularized time evolution changes 
the measure of a volume element in phase space by an 
amount related to (but different from) the variation of 
the Liouville volume. The variation per unit time and 
unit mass of a volume element, measured via /in in the 
sector of phase space containing Nj > particles in flj n 
A, j = 0, 1, . . . , v, can be computed and is, under the 
A-regularized dynamics, 



(l--^r)+f3o(Ko + Uo)) (2.3) 



j>0 



k B Tj{x) v dN. 



as it follows by adding the time derivative of (3q{Kq + Uq) 
to the divergence of Eq.(1.3) (regarded as a first or- 
der equation for the q's and q's) using the expression 
in Eq.(1.5) for Uj. 

Remarks: (1) The dynamics given by the Eq.(1.3) or by 
the same equations with otj = arc different. 
(2) The relation [3 (K + Uo) =(3 {&-± - E j>0 (Uoj - 
Qj)) is useful in studying Onsager reciprocity and Green- 
Kubo formulae, [6]. 
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(3) It is also interesting to consider isokinetic ther- 
mostats: the multipliers ctj are then so defined that Kj 
is an exact constant of motion: calling its value ^Njk R Tj 
the multiplier ctj becomes 



def Qj - Uj 



dN 3 k B T 3 lm 



(2.4) 



with Qj defined as in Eq.(1.5). They have been studied 
heuristically, from the present point of view, in [1]. 

Choose initial data with the distribution fiQ and let 
x — > S^'^x, a = 0,1, be the solution of the A- 
regularized equations of motion with ctj = (a = 0, 
"Hamiltonian thermostats" ) or a;ternatively ctj given by 
Eq.(1.5) (a = 1, "isoenergetic thermostats"), assuming 
t < t A (x). 

Let S^x be the infinite volume dynamics limA^oo 



0,1, provided the limit exists. Let 



x (A.l) {t) d*f s (A,l) x ^ x( A,0) {t) H/^A.O)^ 

x^(t) d Msl 0) x. 



(2.5) 



In the Hamiltonian case the existence of a solution to 
the equations of motion poses a problem only if we wish 
to study the A — > oo limit, i.e. in the case in which 
the thermostats are infinite. For A finite S^ A ^x is well 
defined with /^-probability 1 as in [7]. 

In the thermostatted case the kinetic energy appearing 
in the denominator of ctj, see Eq.(1.5), can be supposed 
to be > with /xo-probability 1. However it can be- 
come at some later time t\(x) (see item (6), p. 2, and 
the example at the end of Sec. III). In the course of the 
analysis it will be proved that with /xo-probability 1 it 

:(A,i) 



18 *a(*) A=^» 

fined. 



oo ; therefore 5 t x is eventually well dc- 



We shall denote (X$ A, ° ; (t),Xp°'(t)) or (S, 



(A,a), 



,(A,a) 



x)i or 



„(A,a] 



(t) the positions and velocities of the particles of 



SI 



(A,a), 



And by Xjf (t) the pairs of positions 

and velocities (gj A ' a ^(i), (^'"^(t)) with q t e Oj. 

Then a particle with coordinates (qi,Qi) at t = in, 
say, the j-th container evolves, see Eq.(1.3), as 

t 



qi(t) =<fc(0) + / qi(t')dt' 
Jo 

q t {t) =e-lo aa ^' )dt 'qi(0) 

Jo m 



(2.6) 



where F t {t) = -d qi (Uj(Xj{t)) + Uj, o(X (t), X, (t))) + 

5j *i(Xo(t)) and Xj(t) denotes X< A ' a) (i) or xj 0) (t). 

The first difficulty with infinite dynamics is to show 
that the number of particles, and their speeds, in a finite 
region A remains finite and bounded only in terms of the 



region diameter r (and of the initial data): for all times 
or, at least, for any prefixed time interval. 

It is convenient to work with dimensionless quantities: 
therefore suitable choices of the units will be made. If 6 
is a prefixed time which is the maximum time that will 
be considered, then 



ipo : (energy scale), r v : (length scale) 
6 : (time scale), v\ - 



^(vel„„) ,27) 



m 

are natural units for measuring energy, length, time, ve- 
locity, respectively. 

It will be necessary to estimate quantitatively the size 
of various kinds of energies of the particles, of a config- 
uration x, which are localized in a region A. Therefore 
introduce, for any region A, the following dimensionless 
quantities: 



(a) N&(x), Nj t &(x) the number of particles of x 
located in A or, respectively, A n % 



(2.8) 



(b) E A {x) d ^ max(^-+V(<z))/Vo 

(c) U A (x) = i _ Vj)/<Po 

(d) Va (x) = max -^J- 

9i£A V\ 



The symbol B(£, R) will denote the ball centered at 
£ and with radius Rr v . With the above notations 
the local dimensionless energy in £>(£, R) will be defined 

def 

asVK(x;£,i?) = E mR) (x)+U B( ^ R) (x) + N B{ ^ R) (x) or, 
more explicitly, 



W(x;£,R) 



def 1 



E 



mqf 

tp(0) ^ \~ 
^ v ' Qi eB(t,R) 



+ (2.9) 



Let log + z d = max{l, log 2 |z|}, g c (z) = (log + z) c and 



def W(x;£,R) 
oc{x) = sup sup 



£ «>9C(^) 



R d 



(2.10) 



If W is the space of the locally finite configurations (i.e. 
containing finitely many particles in any finite region) 
and let Tic C Ti be the configurations with 



(1) £ c (x) < oo, (2) 



A' 



1 

|Anfij| ' 2 2/3^ 



> 



(2.11) 
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for all A = B(0,L) large enough and for Sj,Tj, given 
by Eq.(2.2). Let Nj t \, Uj t \,Kj t \ denote the number of 
particles and their potential or kinetic energy in Ctj n 
A. Each set TLq has /^-probability 1 for £ > 1/d, see 
Appendix A,B. 



III. EQUIVALENCE: ISOENERGETIC VERSUS 
HAMILTONIAN 

It can be expected (and proved here if d — 1,2) adapt- 
ing to the present situation a conjecture proposed in [1], 
that the following property holds for the time evolutions 

configuration x: 

Local dynamics Let d = 1,2,3. Given > 0, with 
Ho -probability 1 then for t G [0, 0], 

(1) The limits x^ a \t) d = lim x (A ^(t) ("thermodyna- 

A^oo 

mic limits") exist for all t < G and a = 0, 1. 

(2) Fort< 0, x^(t) satisfies the second of Eq. (2.11). 

(3) The function t — > x^(t) solves uniquely the Hamil- 
tonian equations in a subspace ofH to which also x^(t) 
belongs (explicit, sufficient, bounds are described in the- 
orem 7). 

Remarks: (a) The limits of ir/ A ' a )(t), as A — > oo, are un- 
derstood in the sense that for any ball A whose bound- 
ary does not contain a particle of x^ (t) the labels of the 
particles of x^ (t) and those of the particles in x( A ' a ^ (t) 
which are in A are the same and for each i the limits 
liia^ QC (ql A ' a \t) , <Z- A '"' ) (t)) exist and are continuous, to- 
gether with their first two derivatives for each i. 

(b) Uniquess in (3) can be given several meanings. The 
simplest is to require uniqueness in the spaces H.^ for 
( > 1/d fixed: and theorem 9 shows that for d = 1,2 one 
could suppose such simpler property. However our result 
is more general and we have left deliberately undeter- 
mined which subspace is meant in (3) so that the deter- 
mination of the subspace has to be considered part of the 
problem of establishing a local dynamics property. The 
generality might become relevant in studying the case 
d = 3, where even in equlibrium there is no proof that 
the evolution of data in £q remains in the same space. 

(c) Recalling the characteristic velocity scale (namely 
vi = ^J2(p(0) /m), the initial speed of a particle located 
in q G M. d , \q\ > r v , is bounded by viy/£gi/ d (q/r v )% ; and 
the distance to the walls of the particle located at q is 
bounded by (y/£ g 1/d (q/r v ) d )- 1/a r v . 

(d) Hence for \q\ large they are, respectively, bounded 
proportionally to [(log \q\/r v )i] i and [(log \q\/r v )i]i : 
this says that locally the particles have, initially, a fi- 
nite density and reasonable energies and velocity distri- 
butions (if measured on boxes of a "logarithmic scale"). 
The theorem 9 in appendix B will show that this property 
remains true for all times, with /xq probability 1. 



(e) An implication is that Eq.(2.6) has a meaning at time 
t = with /io^probability 1 on the choice of the initial 
data x, because £{x) < oo. 

(f) The further property that the thermostats are effi- 
cient: i.e. the work performed by the external non con- 
servative forces is actually absorbed by the thermostats 
in the form of heat Qj, so that the system can eventually 
reach a stationary state, will not be needed because in a 
finite time the external forces can only perform a finite 
work (if the dynamics is local) . 

(g) It has also to be expected that, with /^-probability 
1, the limits in item (2) of Eq.(2.11) should exist and be 

equal to %w- for almost all t > respectively: this is a 

Zjjj 

question left open (as it is not needed for our purposes). 

Assuming the local dynamics property, equivalence, i.e. 
the property x^°\t) = x^(t) for any finite t, can be 
established as in [1]. This is recalled in the next few lines 
of this section. 

In the thermostatted case, with A-regularized motion, 
it is 

= rX^CMX^M (3.1) 
The force between pairs located in f2 , % is bounded by 

def 

F = max \dip(q)\; the numerator of Eq.(3.1) can then be 
bounded by FN y/ Nq y/2Kj/m where A^o is the number 
of particles in Co and A^e bounds the number of thermo- 
stat particles that can be inside the shell of radii Do, 
Do + r v for < t < (having applied Schwartz' inequal- 
ity). _ 

Remark that the bound on A?o exists by the local evo- 
lution hypothesis (see (1) and remark (a)) but, of course, 
is not uniform in the initial data x). 

For < t < and for large enough A, by Schwartz' 
inequality, 

~ y/2K jA (x^)(t)) ~ y/\STj n Wjd/Wj 

having used property (2) of the local dynamics assump- 
tion; letting A — > oo it follows that aj 0. 

Taking the limit A — > oo of Eq.(2.6) at fixed i, this 
means that, with -probability 1, the limit motion as 
A — > oo (with /3j,Xj, j > 0, constant) satisfies 

q t {t) =qi + [ qi(t')dt', q t {t) = ft + / F t (t")dt" (3.3) 
Jo Jo 

i.e. Hamilton's equations, see Eq.(2.6); and the solution 
to such equations is unique with probability 1, (having 
again used assumption (3) of the local dynamics). In 
conclusion 



G 
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Theorem 1: If the dynamics is local in the above sense 
then in the thermodynamic limit the thermostatted evolu- 
tion, within any prefixed time interval [0, 0], becomes the 
Hamiltonian evolution at least on a set of configurations 
which have probability 1 with respect to the initial distri- 
bution jiQ, in spite of the non stationarity of the latter. 

Suppose, in other words, that the initial data are sam- 
pled with the Gibbs distributions for the thermostat par- 
ticles (with given temperatures and densities) and with 
an arbitrary distribution for the finite system in Q with 
density with respect to the Liouville volume (for instance 
with a Gibbs distribution at temperature T and chemi- 
cal potential Ao as in Eq.(1.6)). 

Then, in the thermodynamic limit, the time evolution 
is the same that would be obtained, in the same limit 
A — > oo, via a isoenergetic thermostat acting in each 
container Clj n A and keeping its total energy (in the 
sector with Nj particles) constant and with a density 
equal (asymptotically as A — > oo) to ej. 

The difficulty of proving the locality property (2) can- 
not be underestimated, although it might seem, at first 
sight, "physically obvious" : the danger is that evolution 
implies that the thermostat particles grind to a stop in 
a finite time converting the kinetic energy entirely into 
potential energy. The consequence would be that ctj be- 
comes infinite and the equations even ill defined. As a 
consequence it is natural to expect, as stated in the lo- 
cal dynamics assumption, only a result in /io~probability. 
This can be better appreciated if the following counterex- 
ample, in the Hamiltonian case, is kept in mind. 

Consider an initial configuration in which particles are 
on a square lattice (adapted to the geometry): regard the 
lattice as a set of adjacent tiles with no common points. 
Imagine that the particles at the four corners of each tile 
have velocities of equal magnitude pointing at the center 
of the tile. Suppose that the tiles sides are > r v . If <p(0) 
is large enough all particles come to a stop in the same fi- 
nite time and at that moment all kinetic energy has been 
converted into potential energy: at time all energy is 
kinetic and later all of it is potential. Certainly this ex- 
ample, which concerns a single event that has, therefore, 
probability in /j,q, shows that some refined analysis is 
necessary: the thermostatted evolution a;^' 1 ^*) might 
be not even be defined because the denominator in the 
definition of ctj might become 0. 

It should be stressed that the thermostats models con- 
sidered here preserve even at finite A an important sym- 
metry of nature: time reversal: this ceertainly explains 
the favor that they have received in recent years in the 
simulations. 

Finally a corollary will be that the non dissipative 
Hamiltonian motion and the dissipative thermostatted 
motions, although different at finite volume become iden- 
tical in the thermodynamic limit: neither conserves phase 
space volume (measured with //q) but in both cases the 
entropy production rate coincides with the phase space 
^o^volume contraction. 



IV. FREE THERMOSTATS 

The need for interaction between the particles in order 
to have a physically sound thermostat model has been 
stressed in [8, 9] and provides a measure of the impor- 
tance of the problems met above. 

The above discussion is heuristic because the local dy- 
namics assumption is not proved. However if the model 
is modified by keeping only the interaction ip between the 
test particles and between test particles and thermostat 
particles, but suppressing the interactions between par- 
ticles in the same Qj, j > 0, and, furthermore, replacing 
the wall potentials by an elastic collision rule. I.e. sup- 
posing Uj(Kj) = 0, j > 0, together with the collision 
rule, the analysis can be further pursued and completed. 
This will be referred as the "/ree thermostats" model. 

It can be remarked that in the Hamiltonian case this is 
the classical version of the Hamiltonian thermostat mod- 
els that could be completely treated in quantum mechan- 
ics, [2]. 

Let A„ be the ball B(0, 2") of radius 2 n r v and n> n 



be such that 2 r ' 



> D - 



if N bounds the number of 



particles in the ball Dq + r v up to an arbitrarily prefixed 
time 0, the first inequality Eq.(3.2) and the supposed 
isoenergetic evolution (which in this case is also isoki- 
netic) 



\ aj \<N F< 



N 



< 



N F 



def 



2Kj/m ^Jdk B T 3 lm 



(4.1) 



It follows that, for ( > 1, the speed of the particles 
initially in the shell A„/A„_i with radii 2 n r v , 2 n+1 r v 
will remain within the initial speed by, at most, a fac- 
tor A ±x = e ±ee . The initial speed of the latter particles 
is bounded by, see Eq.(2.10), 



V n = v 1 



n hCd 



(4.2) 



Hence if n(0) is the smallest value of n for which the 
inequality 2 n r ip — V n A < Dq + r v does not hold no 
particle at distance > 2™( e ) +1 r ip can interact with the 
test system. 

This means that 77 < £ c (x) 2 ("( e ) +1 ) d and the dynam- 
ics a;( n > a ) {t) becomes a finitely many particles dynamics 
involving < N + N particles at most. 

From the equations of motion for the Nq + N particles 
we see that their speed will never exceed 



Ve = (V + FN NQ)X 



(4.3) 



if V is the maximum of their initial speeds. In turn this 
means that for n large enough a better bound holds on 

Ctj 



3 ! 



\ajqi 



< 



NqNVqF 



Lo 3 . n 2 dn rdSk B T/m 



^0 



(4.4) 
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with T = min^o Tj and 6 — minj>o Sj and uij. n 2 dn r d 
bounds below (for suitable 0Jj, n ) the volume of flj n A n . 

Hence, for a = 0,1, it is lim^^ (t) = x^(t), 
and also the dynamics is local in the above sense. This 
completes the analysis of free thermostats and proves: 

Theorem 2: Free isoenergetic and Hamiltonian ther- 
mostats are equivalent in the thermodynanic limit 

Notice that essential use has been made of the prop- 
erty that, in absence of interaction among pairs of ther- 
mostat particles and in presence of perfectly elastic walls, 
isokinetic and isoenergetic dynamics coincide: so the de- 
nominators in Eq.(4.1) are constant. 

It would be possible to consider non rigid walls, mod- 
eled by a soft potential ip diverging near them. We do not 
perform the analysis because it is a trivial consequence of 
the analysis that follows. We have chosen the example of 
this section because it pedagogically illustrates well the 
simplest among the ideas of the coming analysis. 



def 

Calling £ = £ 1 / d (x), Eq.(2.10), consider a sequence of 
balls A„ = ,8(0,2"), of radii L n = 2 n r v with n > n : 
so that 2 n °r v > D n + r v and all A„ enclose the test 
system and the particles interacting with it. Given a 
configuration x define N(x; £, R) the number of particles 
in the ball of radius Rr v centered at £ and, given x, 



(1) V n = the maximum velocity in A„ 

(2) p n = the minimum distance to d(Oj n A„) (5 2) 

(3) Af n = max N(x;qi, 1) 

9i6A„ 

After the definition of W, £, the initial speed of a particle 
in the ball A„, n > 1, and its distance to the walls will 
be bounded above and respectively below by v n , Pn with, 
see definitions Eq.(2.7), 



V n = vi (n£)- 



(n£y/ c 



M n < C(n£)i (5.3) 



V. KINEMATICS 

The proof of the local dynamics property will require 
controlling the maximal particles speeds, the number of 
particles interacting with any given one as well as their 
number in any finite region. This will be achieved by 
proving bounds on the local energies W(x; £, R), Eq.(2.9). 

In this section we shall prove bounds at time 0, see 
Eq.(5.3) below as a preparation to the next section where 
we shall use energy conservation to extend the bounds to 
positive time. 

To study general thermostats in dimension d consider 
the dimensionlcss sum of the energy, measured in units of 
tpo, and the particle number in the ball £?(£, R), with cen- 
ter £ and radius r v R as defined in Eq.(2.9) and denoted 
W(x;Z,R). 

The potential ip is superstable so that the number TV 
of points in a region A can be bounded in terms of the 
potential energy U in the same region and of ipo — y(0) 
and Tp > (defined after Eq.(1.4)). This is checked below. 



In fact, 



Eq.(2.9), W > + N ) > 



by the dimensionlcss energy definition in 

17 X ^ ^ 2^Ep N p with thc sum 
running over labels p of disjoint boxes of diameter -re- 
covering A and containing N p > particles (in particu- 
lar: N = J2 p n p)i hence over t < \A\(2Vd/ rip ) d terms. 
By the Schwartz' inequality y/Tp/2p N < \JW I gives a 
bound of the total number of particles in a region A in 
terms of the local dimensionlcss energy W: 



N A < C- 



W 



Tp 



(5.1) 



This is the well known "superstability estimate" (derived 
in our simplifying assumptions of <p > and finite range) . 



under the assumption that the wall potential has range 
r$ and is given by Eq.(1.2); the last inequality is a conse- 
quence of the definition of W and of the above mentioned 
superstability. 

Constant convention: From now on we shall encounter 
various constants that are all computable in terms of the 
data of the problem (geometry, mass, potentials, den- 
sities, temperatures and the (arbitrarily) prefixed time 
6) as in the above Eq.(5.1) which gives a simple exam- 
ple of a computation of a constant. To avoid prolifera- 
tion of labels all constants will be positive and denoted 
C, C, C", . . . ,B,B', . . . or c, c', c", ...,b,b', b", . . .: they 
have to be regarded as functions of the order of appear- 
ance, non decreasing the ones denoted by capital letters 
and non increasing the ones with lower case letters; fur- 
thermore the constants C, . . . , c, . . . may also depend on 
the parameters that we shall name £ or, in Sec. IX, E and 
will be again monotonic non decreasing or non increas- 
ing, respectively, as functions of the order of appearance 
and of £ or E. 

Consider motions, evolving for times < t < 6, or in 
the thermostatted case for < t < min{tA ?i (x), 0}, from 
an initial configuration x following the K n -regularized 
evolution of Sec. 2 with n fixed (see comment (6), p. 2). 
Define 



R n (t) d ^n*+ fv n {s)-, (5.4) 
Jo r V 

where R n (0) — <?i/d(2") = n x l d and V n (s) is the max- 
imum speed that a moving particle can acquire in the 
time interval [0, s] under the A„-regularized evolution. 

The choice of i?„(0) = is made so that it will be 
possible to claim that W(x(0); O, R(0)) < £(x(0))n with 
/Uo-probability 1, see Eq.(2.9),(2.10) and appendix A. 



8 



VI: Energy bound 



The dimensionless quantity R n (t) will also provide a 
convenient upper bound to the maximal distance a mov- 
ing particle can travel during time t, in units of r v , fol- 
lowing the A„-regularized motion. 

It will be necessary to estimate the total energy and 
number of particles in a ball of radius R n (t)r v around £ 
assuming the particles to move with the A„-regularized 
equations. Consider the ball around £ £ M. d of radius 



fl n (t,s) d = f R n {t)+ fY^-ds> 1. (5.5) 

J s r ¥ 

This is a ball whose radius shrinks as s increases between 
and t at speed V n (s): therefore no particle can enter it. 
Abridging x^ n ' a \r) by x(t), this can be used to obtain 
a bound on the size of W(x(t); £, R n {t)) in terms of the 
initial data x(0) = x and of 



W n (x,R) d = f snpW(x n ;^,R). (5.6) 
i 



VI. ENERGY BOUND 



A. Hamiltonian systems 



Considering W(x(s);^,R n (t, s)), for < s < t < 6, it 
follows that 



ds 



W„(x(s);£,i?„(M)) < — V xdl(s),Rn(t,s)) 



q€A„ 



■£( T *T- + 1>M»)) + 5 E v(«W-«'(-)))(6.i) 



«'6A„ 



because the s-derivative of X?(<7( s ): Rn(t, s)) is < since 
no particle can enter the shrinking ball £?(£, R(r, s)) as s 
grows: i.e. X|(?( s )> ^?n(*> s )) cannot increase. 

In the Hamiltonian case a computation of the deriva- 
tive in Eq.(6.1) leads, with the help of the equations of 
motion and setting x^(g(s), R n (t, s)) = X£,q,t,s, to 



if x n denotes the particles of x in A„ = B(0, 2" + 1), i.e 
at distance < r v from A„. 

TTie analysis in the following Sec. VI, VII is taken, with 
a minor adaptation effort, from the version in [10, p. 34] 
of an idea in [11, p. 72] and is repeated here only for com- 
pleteness. 

Let Xz(°:iR) De a smooth function of q — £ that has 
value 1 in the ball R) and decreases radially to reach 
outside the ball £>(£, 2R) with gradient bounded by 
(r v J?) -1 . Let also 



<?eA„ 



(5.7) 



Denoting B an estimate of how many balls of radius 1 are 
needed to cover a ball of radius 3 (a multiple of the radius 
large enough for later use in Eq.(6.3)) in M. d so that every 
pair of points at distance < 1 is inside at least one of the 
covering balls, it follows that W(x;£,2R) < BW{x,R), 
see Eq.(5.5), so that for £ e A„: 



W{x n \Z,R) < W n (x;£,R) < W{x n ;£,2R), 
W n (x;Z,R) < BW n {x,R) 



(5.8) 



Although W has a direct physical interpretation W turns 
out to be mathematically more convenient and, for our 
purposes, equivalent by Eq.(5.8). 



—W n {x(s);^R n (t,s)) < ?(*)*(3(*))X£,«,M 



qei2 

q{s)d q ip(q(s) - q'(s)) 



q'eA„ 



2(fi 



(6-2) 

where the dot indicates a s-derivative and it has been 
kept in mind that positions and velocities of the par- 
ticles outside A„ are considered, in the A„-regularized 
dynamics, to be time independent. 

Since the non zero terms have \q(s) — q'(s)\ < r v and 
the derivatives of \ are < (r v R n (t, s))^ 1 and \q\, \q'\ < 
V n (s) = r v \R n (t,s)\ it follows, setting F = max.(\dip\ + 

i*n, 



V(^);(,fi n (M))<— ^(xW;C,-) 

ds (f r v 



Fr v \R n {t,s) 

ip R n (t,s) 



■BW n (x(s);S,2R n (t,s) + l) (6.3) 



< B 2 Fvi^l^it^ + 1)Wn{x{s);Rn{ti s)) 

(flO Vl Rn(t,S) 



where W n (x; R) is defined in analogy with Eq.(5.6). 
Eq.(6.3), R n {t,s)/R(t,0) < 2 and W n (x(s), £, Rn(t, a)) 
< W n (x(s),R n (t,s)) imply the inequality 



W n (x(s),R n (t,s)) < e r W n (x(0),R n (t,0)), (6.4) 



with r =' log 2 + 5^1) e b. 



VII: Entropy bound and phase space contraction 



9 



B. Thermostatted systems 



In the thermostatted dynamics case, Eq.(6.2),(6.3) 
have to be modified by adding to the r.h.s the work, per 
unit time, done by the thermostatic forces (measured in 
units of po) which is, see Eq.(1.5), 



£,*g«(*)*Us(s)) 



g»eB(|,fl„(t,s)) 

(6.5) 

here the * means restriction of the sum to the particles 
Qiit) in Clj and within distance r v from the boundary 
of do, which is a ball of radius D . Then Eq.(6.5) is 
bounded by 



N — W n {x{s),R n {t,s)) 
fa 



(6.6) 



because the bounds of the various terms in Eq.(6.5) can 
be obtained as: 

(a) the ^2 i <ii(s)F i (q(s)) by (Schwarz' inequality) < 

NoF(52* <7?) 3 (£i I) 5 with F = max \dip\. Leaving aside 
the factor FNq the rest is bounded above proportionally 
to W n (x(s),R n (t,s)). 

(b) the kinetic energy in the last sum in Eq.(6.5) is also 
bounded by the total kinetic energy and is compensated 
by the denominator. 

(c) No is bounded in terms of £ (and constant in time) . 

Therefore, given n, the new inequality which replaces 
Eq.(6.3) in the thermostatted case gives a bound of the s- 

derivative W n d = j^W n (x(s);£, R n (t, s)) in terms of, see 

Eq.(5.6),(5.8), W n = f sup c W n (x(s); £, R n (t, a)), namely 



W n <C(^+(No + l)^)w n 

R^lt.s) G3n / 



Fvi 



Rn(t,s) Ifo J 

The second addend in Eq.(6.7) is bounded, by (c) 
above. The differential inequality Eq. (6.7) (together with 
Eq.(5.6),(5.8)) implies, for suitable C , C 2 functions off: 

W n (x(t),R n (t)) < C'W n (x(0),R n (t)) < C 2 R n (t) d 

(6.8) 

for t < mm{t\ n (x),Q} (see p. 2, item (6)). Hence the 
bound on the speed, for instance for d = 2, Vn ^ < CRn 2 
with R n = R„{e). 



C. The bound 

Therefore the following energy bound holds: 

Theorem 3: For a suitable constant C the follow- 
ing energy bound holds for the A„ -regularized Hamil- 
tonian or thermostatted dynamics and for t < 9 or 
t < min{*A„(^), ©}, respectively, 



W n (x^°\t),R n (t)) < C 2 R n {t) d (6.9) 

and n large enough and C > (depending only on £). 

Remark: the inequality holds for all d's and in the Hamil- 
tonian case the constant C can be taken £ -independent, 
see Eq.(6.4). 

By the definition Eq.(2.9) of W it follows that 



V n (s)<v 1 CR n (s) i - 



(6.10) 



and going back to Eq.(5.5), and solving it, R n {t) is 
bounded proportionally to i?„(0) = n 1/<d . 

Calling p n (t), V n (t),Af n (t) the quantities in Eq.(5.2) for 
S[ n,1 ^x the following bounds can be formulated simulta- 
neously for the thermostatted and Hamiltonian cases and 
extend to positive times the estimates in Eq.(5.3). 

Theorem 4: If d < 2 the maximal velocity V n (t) and 
the maximal displacement R n (t) r v , in the K n -regularized 
motion up to time t < min{6, t^ n (x)} in the thermostat- 
ted case or t < in the Hamiltonian case satisfy: 



R n (t)<Cni, V n (t) < v x C n 1/2 

P n(t) ^C- 1 ^^, N n {t)<Cn^ 



(6.11) 



with C > (a suitable function of £). 

The use of Eq.(6.10) to solve Eq.(5.4), which then im- 
plies Eq.(6.11), will force the restriction on the dimension 
to be d < 2. However if the Eq.(6.11) are assumed also 
for d = 3 the following analysis and results would hold 
true also for d = 3, except a "minor" modification of the 
extra result in theorem 9 in appendix C; therefore we 
keep in the following the dimension d generic. 



VII. ENTROPY BOUND AND PHASE SPACE 
CONTRACTION 

The kinetic energy density Kj iTl (x)/\Qj fl A n |, i.e. the 
kinetic energy contained in the j-th thermostat at time t 
in the A„-regularized thermostatted or Hamiltonian dy- 
namics divided by its volume, will initially have a value 
as close as wished to Sj^ksTj for n large enough. So 

that K 



~ Kj 2 dn if n is large enough: because this is 
an event which has probability 1 with respect to jj,o- 

Therefore if the infinite volume ^o-average value of 
Kj - n ^ is denoted Kj2 nd then there is n(x) such that 



^^l>l K 2 nd , Vn>n(x). (7.1) 
fo 2 

where k = min^oKj. The notation x(t), x(s), . . . will 
be temporarily used below for simplicity instead of 
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x^ n ' a \t), x( n ' a )(s), .... The equations of motion arc now 
Eq.(1.3), or (2.6), with a = 1. 

With /^-probability 1 it is £(x),n(x) < oo. Let Ee,h 
be the set of configurations for which £(x) < E, n(x) = h. 
Let Vh(x) be the smallest n > h such that the event 
min j>0 Ki '^* m = \k 2 nd is realized for t = t nth (x) < 9 
but not earlier. 

Remark: the definition of Ee,Ii implies that the time 
t\ n (x) (sec comment (6), p. 2) is, for x G H^f,, certainly 

> tn,h{x)- 

Let ^E,h,n be the set of the x G z-E,h, Vh(x) = n. It 
will be shown that Y^ n >h W>(^E,h,n) < °° : hence with 
/^■probability 1 a point x will be out of ^E,h,n for all 
n large enough and Eq.(7.1) will hold for alH < O with 
/^■probability 1. 

For x G ^E,h,n, consider Qj as in Eq.(1.5). Notice that 
Qj has the physical interpretation of the heat ceded per 
unit time to the j-th thermostat by the system. 

If TVo is the number of particles in Qq, A/" n (i), V n (t) are 
as in Eq.(6.11), F is the maximum of \d<p\, and N is the 
maximum number of particles within r v of the bound- 
ary of Q , then \Qj\ < N JfFV n (t) < C'ni+i because 
V n (t) is bounded by Eq.(6.11) proportionally to n 1 / 2 , N 
is bounded proportionally to the (finite) number of balls 
of radius r v needed to cover the ball of radius Dq + r v 
times the N n (t) in Eq.(6.11) (i.e also proportionally to 
rja): the estimates hold for t < 8 in the Hamiltonian 
case and for t < th, n (x) < 6 in the thermostatted case.. 

The phase space contraction a = a(x), see Eq.(2.3), is 



a = J2 dN i 

j>0 



3 2K jtn (x,ty dN 3 



ief 



) + /? Qo (7.2) 



in the thermostatted case, if Qq "= — J2 q en ^j>o ^-o' 
dx Uoj{X-o, Xj) + X • ^; in the Hamiltonian case, it is: 



(7.3) 



3>0 



It has to be kept in mind that there is contraction of 
phase space even in the fully Hamiltonian case when Li- 
ouvillc's theorem holds (already for the regularized dy- 
namics): this is no contradiction because the phase space 
volume is measured with the volume defined by pq. 
Since Qo can be estimated in the same way as Qj, 

above, it follows that the integral | f* h ' n{x) a(x^ n ' a \t))\ 
is, in either cases, also uniformly bounded (in E E ^ h ) by 



(707J, 



Vn > h 



(7.4) 



where a is a suitable non decreasing function of E. This 
is a first entropy estimate; it is rather far from optimal 
and it will appear, Sec. IX, that it will be essential to 
impove it. 



Therefore a volume element in Es,h,n contracts at most 
by e _cre ™ on the trajectory of /io-almost all points x G 
^■E,h,n up to the stopping time th, n {%)- 

Effectively this means that the distribution p can be 
treated as an invariant one for the purpose of estimating 
the probability that the kinetic energy, in Qj n A„ of the 
initial data x G ^E,h,ni in the time th, n (x) grinds down 
to § 2 dn (i.e. to half (say) of the value K2 nd to which it 
is intially very close, if n is large). The estimate can be 
carried out via the technique introduced by Sinai, [12], 
which has been applied in [11, 13]; for completeness sec 
the following appendix D. 

Let D = D n be the set of the x G ^E,h,n which satisfy 
Kj t \ n (x) = ^n2 dn for a given j > while Kj>,\ n (x) = 
\n2 dn for / > OJ'^j. 

Recalling the DLR-equations, [14], and the classical 
superstability estimate on the existence of b > such 
that p n < e~ b2 " bounds the probability of finding more 
than p2 dn particles in A„ n Qj if p is large enough (e.g. 
p > maxj Sj), the probability po{^E,h,n) can be bounded 
by p n (summable in n) plus 



po(dq'dq') 



P 2 na 



-(0jU An , j (q,q')-X j I) dq 



1=1 



■ e- f> > p 'P ld - 1 Pu;(ld) 



(7.5) 



±k2*\ 



where q = (qi,...,q t ) G (fi^- n A n )', P 2 
UA n ,j(q,q') is the sum of ip(q — q') over the pairs of 
points qi,q' e G Qj fl A„ plus the sum over the pairs with 
qi G A„ n flj , q' t £ A„ n Qj , and 

(1) Z\ n ,j(q') is the partition function for the region A n n 
Qj (defined as in Eq.(2.1) with the integral over the q's 
extended to A„ n Qj and with the energies U\ n (x, z))\ 

(2) The volume element P ld ~ x dP has been changed to 
p id-ip dT = p id-2p dT where p ig a short hmd for 

E 9 , 9 ' ;(?e A„ \9MQ - «') I + E, g a„ 13^(0)1 so that PP is 
a bound on the time derivative of the total kinetic en- 
ergy P 2 contained in A„ evaluated on the points of D n 
(the latter is 2PP = \ J2i,j m e\ n dt f(<l* - Qj)(.Qi - Qj) + 
E 9£A „ d q Mq)q\ hence < P 2 P). 

(3) uj(ld) is the surface of the unit ball in M. ld . 

(4) The factor e a6n takes into account the entropy esti- 
mate, i.e. the estimate Eq.(7.4) of the non-invariancc of 
Mo- 

The integral over r in Eq.(7.5) gives a factor and 
the integral can be trivially imagined averaged over an 
auxiliary parameter e G [0,e] with e > arbitrary (but 
to be suitably chosen shortly) on which it does not de- 
pend at first. Then if P is replaced by (1 — e)P in the 



the 



exponential while P ld 1 is replaced by — — ^ p2 „ d _ 1 

average over e becomes an upper bound. Changing e to 
Pe (i.e. hence de to = the bound becomes the 

^o-average 
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\n2 nd 



2e^ p2 " d ~ _ 2e^ 2nd - 2 i i 

eK2 d ™ ( PXk '^»° = e K 2 dn ^ ^° ' (7.6) 

where x £jK is the characteristic function of the set {(1 — 
e) 2 § 2 dn < Kj < § 2 dn }. The inequality is obtained by a 
bound on the first average, via a superstability estimate, 
proportional to 2 2dn and by the remark that the second 
average is over a range in which K shows a large devia- 
tion from its average (by a factor 2) hence it is bounded 
above by e~ fc2 ™ with b depending on k but independent 
on £ for n large. Therefore fixing e small enough (as a 
function of k) the bound holds with suitable B,b > and 
is summable in n (and of course on j > 0). 

Hence, fixed h, with /ito-probability 1 it is Kj >n > 
(by Borel-Cantelli's theorem) for all n large 
enough and j > 0. As mentioned after Eq.(7.1) this 
means that for all i < it is Kj tTl > ^n2 nd for all n large 
enough, with /i probability 1. Therefore the bounds in 
Eq.(6.11) can be assumed, with /j -probability 1 also for 
the thermostatted and Hamiltonian dynamics and for all 
n large enough. 

Theorem 5: With /iq -probability 1 the phase space 
contraction cr(x) admits a bound a(x) < an for all 
times t < 6 for the A n -regularized thermostatted or 
Hamiltonian dynamics. Furthermore the kinetic energy 
Kj,n(x i - n ' a ) (t)), a = 0, 1, in the j-th thermostat remains 
> n2 nd for all n > n(x) for suitable K,n(x) > 0. The 
constants a(x), n{x) depend on x only through £(x). 

This proves item (2) of the local dynamics property. 



VIII. INFINITE VOLUME HAMILTONIAN 
DYNAMICS 



The existence of the dynamics in the Hamiltonian case, 
a = 0, will be discussed first, proving 

Theorem 6: If x £ Ti-i/d the thermodynamic limit evo- 
lution x(°\t)i — linin^oo x( n '°\t)i exists. 

The following proof reproduces the proof of theorem 
2.1 in [10, p. 32] for d — 2, which applies essentially unal- 
tered. Define 



5^n) d M\at fi) (t)-at +lfi) (t)l 

Uk(t, n) d = max Si(t, n), 
<J;eA fc 

then, for a = 0, Eq.(8.1) yields 



Si(t,n) < I -dr {F^6i(r. 



(8.2) 



F'2(5 J -(r,n) + 5 i (r,n))} 



(8.3) 



where the sum is over the number Af n of the parti- 
cles qj{r) that can interact with qi(r) at time r; F' — 
max \d 2 (p\ is the maximal gradient of the intcrparticle 
force; F' = Ca(a + l)^n^3^ bounds the maximum 

gradient of the walls plus the stirring forces, the bound 
follows from Eq.(6.11). The number N n is bounded by 
theorem 4, Eq.(6.11), by N n < Cn 1 / 2 for both x^(t) 
andx( n+1 ' )(r). Let 



/=(1 + ^), 2 k ^^2 k + r n (8.4) 

where r„ is the maximum distance a particle can travel 
in time < 6, bounded by Eq.(6.11) by Cr^n 1 / 2 . Then 



It remains to check that also the n — > oo limit of 
the dynamics exists in the sense of the local dynam- 
ics assumption (i.e. the existence of the limit x(t) = 
x (o) ^ d ^f ]] mn _^ oo x (nfi) (^ an( j a sm table form of its 
uniqueness) . 

The equation of motion, for a particle in the j-th con- 
tainer (say), can be written both in the Hamiltonian and 
in the thermostatted cases as 



aaj(x ( "- a) (s))ds 



q ^ a) (t)= qi (0)+ [\e-Io T - 
Jo 

+ (t-T)e~ T aa > {x(n ' a) {s))ds ft (x^ ( T ))) 



«i(0) 



(8.1) 



dr 



where the label j on the coordinates (indicating the con- 
tainer) is omitted and fi is the force acting on the selected 
particle divided by its mass (for j = it includes the 
stirring force). The Hamiltonian case is simply obtained 
setting a = while the thermostatted case corresponds 
to a = 1. 



u k (t,n) <Cn n r u kl (s,n) ds 



Jo r v 8 



(C is a function of £ as agreed in Sec.V). Eq.(8.5) can 

be iterated £ times if 2 k + C in 1 1 2 < 2", i.e. 1= 

which is I > c2 n l 2 S k<n for n large. 

By Eq.(6.11) u k (t,n) is < Cn 1 / 2 so that for n > k, 



Uk{n,t) < c , (ni) 



V2< C2 - 2 " /c (8.6) 



for suitable C , C, c > (n-independent functions of £). 
Hence the evolutions locally (i.e. inside the ball A*,) be- 
come closer and closer as the regularization is removed 
(i.e. asm oo) and very fast so. 
If qt (0) G A k , for n > k it is 



n—k 
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showing the existence of the dynamics in the thermody- 
namic limit because also the inequality, for n > k, 



/j(n.O) 



(*) 



1 Wi<C2- 2 '° (8.8) 



follows from Eq.(8.6) and from q (n ^(t) - g("+ 1 >°)(t) = 
Jo (fi(Q {nfi) (r))dr - fM {n+m {r)))dr. Or, for n > k, 



>.0)/ 



-c2 



Lti/2 



(8.9) 



calling \xi - x[\ = fa- q^/vi + fa - q'i\/r v . 

Hence the proof of the existence of the dynamics in 
the Hamiltonian case and in the thermodynamics limit 
is complete and it yields concrete bounds as well, i.e 

Theorem 7: There are C(£),c(£) _1 , increasing func- 
tions of £ , such that the Hamiltonian evolution satisfies 
the local dynamics property and if qi (0) € 



\q < r'° ) (t)\<v 1 C(£)ki, 

>,0) 



distance^ ■>{t),d(U j n j n A)) > c{£)k~^r^ 
Ni{t,n) < C(£)k^ 2 
\xt fi \t)-xf\t)\<C{£)r v e-<^ 



(8.10) 



nd/2 



for all n > k. The x^(t) is the unique solution of the 
Hamilton equations satisfying the first three of Eq.(8.10). 

The uniqueness follows from Eq.(8.3) and we skip the 
details, [15]. 

It would also be possible to show the stronger result 
that x^\t) e Ti-i/d'- but for the proof of theorem 1 
the theorems 5,6,7 are sufficient, hence the proof of the 
stronger property is relegated to theorem 9 in the ap- 
pendix. 

The corresponding proof for the thermostatted evo- 
lution will be somewhat more delicate: and it will be 
weaker as it will not hold under the only assumption that 
£{x) < oo but it will be necessary to restrict further the 
initial data to a subset of the phase space (which however 
will still have /io-probability 1). 

Remark: An immediate consequence is that the entropy 
production cr(x), see Eq.(7.3), is estimated by a constant 
s(£) in the A„-dynamics for n large: i.e a much better 
estimate than the growth bounded by a power of n, see 
Eq.(7.4), implied by Eq.(6.11). 



IX. INFINITE VOLUME THERMOSTATTED 
DYNAMICS 

Eq.(8.1) will be used to compare the Hamiltonian and 
the thermostatted evolutions in Clj with the same initial 
data assuming that the initial data satisfy theorem 5, 
Sec. VII. We shall see that the problem will reduce to 



obtain a better estimate of the entropy production, i.e. 
better than proportional to n, as in Eq.(7.4). 

Fixing once and for all n > to be smaller than the 
minimum of the kinetic energy densities of the initial x in 
the various thermostats (which is x-independent with /j,q- 
probability 1), the problem can be solved by restricting 
attention to a suitable subset of the set Xe C H-i/d- 

X E = f {x\£(x) <E; min ^ > n2 nd \ (9.1) 

In this section (and in the corresponding appendix E) 
the constants C, C" , . . . , c, c', . • ■ will be functions of E as 
stated in Sec.V. 

Consider the bands of points £ at distance pn (£,) from 
the boundary <9fi of £1 within r v or 2r v 



K d = {q:pnM < 



A* 



ief 



{<? = Pn„(q) < 2?>} 



(9.2) 



By the result in Sec. VIII there are M and V (which de- 
pend on E) so that for all x € Xe and with the notations 
Eq.(2.8), for n large enough: 



maxJV A (Si n ' 0) x) < M 
t<e 

max V A . (S { t n ' 0) x) <V-1 
to 



(9.3) 



Define for x the stopping times 



T M ,V;n(x) d = max{t < mm{t An (x), 0} : Vt < t, 

Na^S^x) < M, VaAS^x) < V) (9.4) 

Let the cube with side r v centered at a point £ in the 
lattice r v Z d , and using the definitions in Eq.(2.8), let 



I, H def ma^(N C( (x), v /E Ci {x)) 

\\x \\n — max — — — . ly-OJ 

with i < A < 1. Split X E = AUB where 



A= f {xeX E : max \\Si n ''>x\\ n < (logn) A } 

t<T M ,V;n(x) 



?(«>!) „ 



B= f {x e X E : max WS^'x^ > (logn) A }. 

t<T M ,V; n(x) 

(9.6) 

Fixed, once and for all, 7 > arbitrarily 

Theorem 8: In d = 1,2 there are positive constants 
C,C',c depending only on E such that for all n large 
enough: 

(1) if x & A then T M y-,n{x) = 9, S\ nfi) x and s\ n ' 1] x 
are close in the sense that for qi(0) <E A(i ogn ) 7 
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k KD (t) _ g Ko )(t)| < (7rye -(iogr t rc 5 

|«i n,1) (t)-«| n,0) (t)| < C^e"^")^. 
(%) i/ie sei B has ^-probability bounded by 



{B)< Ce~ c( - losn)2X+c ' M2y . 



(9.7) 



(9.8) 



Remark: Since A > 1/2 Eq.(9.8) will imply (by Borel- 
Cantelli's lemma) that, with ^-probability 1, eventually 
x is in A and therefore, as soon as Eq.(9.7) will have 
been proved, the thermodynamic limits of qr(">°)(t) will 
coincide for t < 9, a = 0,1: concluding the proof of 
theorem 1 as well. 



A. Check of Eq. (9.8) 

We begin by defining the surface, see it symbolic de- 
scription in Fig. 2 below: 



E T = f {xeE E \\\S^x\\ n >(\ogn)\ 
Vt<r, ||5 t ( "' 1) x||„<(logn) A 

T M ,V;n(x) > t} 



(9.9) 



Recalling the definition of the existence time t\ n (x) in 
Sec. I, p. 2, item (6), we remark that for x G He, t\ n (x) > 
so that Si n,1 ^x, r < is well defined. 

Moreover S , |"' 1 ' ) S T is contained in the surface E' (piece- 
wise smooth) of points x for which S^"' 1 ^ is well defined 
for t < near and ||5j n,1 ^x|| crosses from below the 
value (logn) A at time t = 0. 

With the above geometric considerations (see Fig. 2) 



the set E d = U T < e Si"; iJ a; G E T C E' and for x G E we 
define 9(x) = max T <e {t | S^^x G E T }. 




Fig. 2: TTie horizontal "line" represents E' ; </ie "curve" rep- 
resents the points S^q^x, i.e the set of points in Xe which 
in time 9(x) reach E' and determine on it a subset E; the in- 
complete ( "dashed") lines represent the "levels" E T ,E T /; the 
missing parts are made of points which are not in Xe but have 
an "ancestor" in Xe; the vertical line represents the trajectory 
of the point S^g^x. 

The E T are represented in Fig. 2 as dashed lines to re- 
mind that it might be that the trajectories that reach the 



surface E' will have a value of £ (S^' 1 ^ x) > E or a kinetic 
energy < K2 nd /ip , see Eq.(9.1), at some t G (0, 0(x)). 

Remark: We can also say that E is the subset of the 
surface E' consisting of the points of E' that are reached 
by trajectories of points y G S# within a time interval 

< T M ,V;n(y)- 

In the notations of Appendix D, E is the base and 6(x) 
the ceiling function. We then have 



Mo(#) < / fJ-o,n'(dy) I dtw{y)e 
Js' Jo 



0(y) 



c{y,t) 



(9.10) 



where /l«o,e' denotes the projection of (Xq on E', w = \v x ■ 

n x \ and a(y,t) d = J° t a(y^ n ^(t'))dt' is the phase space 
contraction. By the definition of E T it then follows that 
T M ,v-,n{x) > 6{x), Vie S and \a(y,t)\ < CM 2 V. 

Writing k^ for the smallest integer > (log n) A t/A (£/?>), 
^o,S' almost surely, E' splits into an union over £ G A„ n 
r v Z d of the union of <S>| U <S| , where 



S\ ={y G E' : \y n Q| = k 6 , \y n 5C ? | = 1} 
5| ={?y G E' : y n Q 9 («, g), g) = £J 



(9.11) 



HE, 



def 



((logn)V(£Av)r 



Both yto.s'(5|) and /ito.E'G^f) are bounded by 



Mo 



E'GS) < Ce CM2y V^e- c [( log ") A ^(«/ r ^l 2 , (9.12) 



(for suitable C, c, functions of E). The proof of Eq.(9.12) 
does not involve dynamics but only classical equilibrium 
estimates, the details are expounded in Appendix E. Sum- 
ming (as A > 1) over £ G fij H A„ the Eq.(9.8) follows. 



B. Two remarks 

To prove item (1) and Eq.(9.7), thus completing the 
proof of theorem 8, we shall compare the evolutions 
x^ n ' a \t) with a = 0, 1, same initial datum x G A and 
t < Tm,V;ti(x), the latter being the stopping time defined 
in Eq.(9.4). We start by proving that there is C > so 
that for all n large enough the following holds. 

Lemma 1: Let 7 > 0. For x G A and h > (logn) 7 ; then 



\qt' 1 \t)\<Cv 1 (h logn) A , 
\qt A \t)\ <r v (2 h + C(h logn) A ). 

for qi(0) G A h and t < 0. 

Proof: if x G A and t < TM,v-, n {x) then 



(9.13) 



4 n,1 \t)\ < «i((logn)log + 



V2r 



(9.14) 
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implying: \q^' ' {t)\ < r(t)r v if r(t)r v is an upper bound 
to a solution of Eq.(9.14) with = replacing < and ini- 
tial datum ^"'^(O)! < 2 h r v . And r(t) can be taken 
r{t) d = 2 h + 2 Vl ((logn) log + 2 h ) X ±, for t < 6, provided 

((logn)log + r(e) + V2) A < 2((logn) log + 2 h f (9.15) 

which is verified for all n large enough, because ^ hio ^ h n ^ 
vanishes as n diverges (keeping in mind that h > 

(logn) 7 ). Thus |«| n,1) (i)| < fyr(t), hence |«j n,1) (t)| < 
r v Cr(t) for alU < T M ,v- n (x), i.e. when Eq.(9.14) holds: 
and the lemma is proved. Fix 7 > 0. 

Lemma 2: Let TV and p be the maximal number of parti- 
cles which at any given time < O interact with a particle 
initially in A^ +1 and, respectively, the minimal distance 
of any such particle from the walls in either dynamics 
and for t < Tm,v-, n ( x ) ■ Then 

Af< C(fclogn)\ p> c{k\ogn)- 2X ' a (9.16) 
for all integers k > (logn) 7 . 

Proof: As a consequence of lemma 1 and of theorem 7 
the following properties hold for all n large enough and 
all t < TM,v-,n{x), both for the Hamiltonian and the 
thermostatted evolutions. 

(i) for all Qi e Afe + 2 and a = 0, 1, 

max \q\ n ' a) {t) - m\ <Cr v {k logn) \ (9.17) 

t<T M ,V;n(x) 

(ii) particles € Afe do not interact with those ^ A^+2; 

By Eq.(9.17) we see that if q, e A fe+ i then q- n ' a) {t) e 
Ak+2 so that, by the definition of the set A and by 
Eq.(9.6) (or by theorem 7 in the Hamiltonian case, re- 
calling that A > 1/2), Eq.(9.16) follows. 

C. Check of Eq.(9.7) and comparison of 
Hamiltonian versus thermostatted motions 

def 

We have now all the ingredients to bound 5i(t,n) = 

— 1i n '°\t)\- Let fi be the acceleration of the 
particle i due to the other particles and to the walls. By 

Eq.(9.16) if ft e Afc+i, \fi\ <C(klogn) r >',r } '= f 2d\(l + 
i) so that, subtracting the Eq.(8.1) for the two evolu- 
tions, it follows that for any e Afe+i (possibly close to 
the origin hence very far from the boundary of Afc if n is 
large, because k > (logn) 7 ) 

Si(t,n) < C(fclogn)"V Ild (9.18) 

+ e /Vi(9 (n,1) «)-/i(«/ (n,0) W)|dr. 
Jo 



because, recalling Eq.(9.1), |aj| is bounded proportion- 
ally to 2~ nd . 

Let £ be a non-negative integer, ki such that 

2 k e =2 k +lC(klogn) x (9.19) 

(see Eq.(9.17)) and Uk e (t,n) the max of 5i(t,n) over 
< 2 fc «. Then by Eq.(9.18) and Eq.(9.16) and writ- 
ing n" = / 2dA(l + |), 

Uk ^ n) <C(k log n) r '' 2-" d (9.20) 
r v 

+ C(fclogn) ?) / ± — . 

for £ < £* = 2 k I ((k log n) A C), the latter being the largest 
£ such that 2 fe « < 2 fc+1 . By Eq.(9.20) and Eq.(9.16) 

u k (t,n) < e c{k losn)v " C(klognf2- dn (9.21) 

Thus u k (t,n) is bounded by the r.h.s. of the first of 
Eq.(9.7); analogous argument shows that also the veloc- 
ity differences are bounded as in Eq.(9.7) which is thus 
proved for all t < TM,v-, n {x). 

On the other hand given qi(0) with \Qi(0)\/r v < 2 k ° 

it is, for n > e k a' 1 and i fixed, (^"'^(i) — < 
U(io g „)7 (t, n) < C2 dn l 2 , i.e. for n large g- n,1 ^(i) is closer 

than r v to q\ n '°\t). Hence, remarking that we know 
"everything" about the Hamiltonian motion we can use 
such knowledge by applying Eq.(9.7) to particles which 
are initially within a distance r v 2 k ° of the origin, with 
ko fixed arbitrarily, 

Therefore the number of particles in qf"' 1 ^ (t) which arc 

in A* is smaller than the number of particles of g-™' ' (t) 
in A** which is bounded by M. An analogous argument 
for the velocities allows to conclude that Eq.(9.3) hold in 
A*, also for the thermostatted motion (a = 1, being valid 
for the Hamiltonian motion in the smaller A*, given the 
closenes of the positions and speeds), Tuy-.n = 6 with 
/io-probability 1. 

Applying again Eq.(9.7) the proof of theorem 8 is com- 
plete: with 7 = 2 (but any 7 > would also lead to a 
corresponding result). 

X. CONCLUDING REMARKS 

Equivalence between different thermostats is widely 
studied in the literature and the basic ideas, extended 
here, were laid down in [16]. A clear understanding of 
the problem was already set up in comparing isokinetic, 
isoenergetic and Nose- Hoover bulk thermostats in [16], 
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where a history of the earlier results is presented as well, 
see also [1, 17]. 

There are, since a long time, studies of systems with 
free thermostats, starting with [2]. Such thermostats are 
somewhat pathological and may not always lead to the 
stationary states that would be expected: as exemplified 
in the case of simple spin chains, [18, 19]. More recently 
similar or identical thermostat models built with free sys- 
tems have been considered starting with [20]. 

The case of dimension 3 is very similar: it is not diffi- 
cult to prove, that the key bounds (6.9) hold; however a 
naive application of the ideas developed in [10] to prove 
that R n (t) satisfies Eq.(6.11) is not possible. 

Isokinetic thermostats should be treated in a very sim- 
ilar way, [1]: the extra difficulty is that the entropy pro- 
duction in a finite time interval receives a contribution 
also from the time derivative of the total energy of the 
reservoirs, [1], and further work seems needed. 

More general cases, like Lennard-Jones interparticle 
potentials are more difficult, see [21]. Finally here the 
interaction potential has been assumed smooth: singular- 
ities like hard core could be also considered at a heuristic 
level. It seems that in presence of hard cores plus smooth 
repulsive potentials all estimates of Sec. V, VI are still 
valid but the existence of the limiting motion as A — > oo 
remains a difficult point because of the discontinuities in 
the velocities due to collisions. 



XI. APPENDICES 
A. Appendix: Sets of full measures 

There are cq and Rq and a strictly positive, non de- 
creasing function 7(c), c > c , so that Vc > c ,VR > Ro, 

p (W(x, 0, R) > Ci? d ) < e -7(")* d (ll.l) 
If g : Jj d — > R + , g(i) > 1, c > c , the probability 

^. eZ ,^ g{i) W(x;i,r)<cr d ) (11.2) 
is > 1 — VI rpd „ ... e -7 ^ 1 * with the sum being 

— ^te£ ,r>g(i) ° 

bounded proportionally to the sum VJ ie ^ £i e~ 7 ^ 9 ^ 

which converges if g(i) > c'(log + |i|) 1 ' d , with c' large 
enough. 

B. Appendix: Choice of R n (t) 

The proof of the inequalities Eq.(6.4),(6.9) yields Vt < 
9 that W(S { t n) x, R) < cW(x, R + f* V n ) provided R is 
R+ f * V„ 

such that — ^ — < 2, which is implied by R > R + 
Jq V n (s)ds, Rq > 0. The maximal speed V(t) at time t is 



bounded by V(t) < v 1 ^2W n (S ( t n, x,R). Choosing R = 

R n (0) = n3 we get V n (t) < C"uii?„(t) d / 2 < Cv x n^: such 
choice is the weakest that still insures that the set of 
initial data has W(x, 0, R)/R d finite with /i -probability 
1, see appendix A. 



C. Appendix: The Hamiltonian motion is a flow in 

C > § 

The following theorem is obtatined by a straightfor- 
ward adaptation to the case d = 1,2 of theorem 2.2 in 
[15]. 

Theorem 9: Let d = 1,2, ( > 1/d, E > 0. Then, given 
any 9 there is E' (depending on Q,Q,E) so that for all 
x such that £q{x) < E 



£d s t x ) ^ E '> f° r aiit<e (11.3) 

so that the evolution x — > S^x is a flow in all spaces 

n c , c > i/d. 

So far, for the sake of definitcness, ( = 1/d has been 
assumed: therefore in the following proof the quantity £ 
has to be intended equal to 1/d; however £ is left unde- 
termined because the proof would still hold for arbitrary 
C > if larger values were consistently assigned to 

it since the beginning of this paper, under the essential 
restriction d < 2. 

Proof. Let x[ n ' ^ = S^ n,0 ^x and consider 

W(x { ?\t;,p), for p>(log + (|CHAv))< (H-4) 

with W defined as in Eq.(5.7) with no restriction in the 
sums over q,q'. Let — 1 be the smallest integer such 
that A„ 4 _! contains the ball of center £ and radius pr v . 
Then Vt < 6 



W(xi°U,p)<W n M ni ' 0) ^P) (U5) 
+ \W(x[ \z, P )~W ni (x ( ^ 0) ,t,p)\ 

The motions x|™ ? ' ' and x^ are very close for all 
points which initially arc in A„ s _i: by Eq.(8.6)-(8.8) 
the difference of positions and velocities are bounded by 
Cexp-c2™«/ 2 . 

Setting X?(?) equal to the smoothed characteristic 
function Xc(<Zj P) introduced in Eq.(5.7), X?(9i™'°^ (*)) an< i 
Xiili (t)) force their arguments to be within if 
2 fei = 2 "«-i + p < 2"«. Hence the inequality 
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•("5.0), 



\W(x\ 0) ,Z,p)-W n( (x. 

< ce-™ «E^k K ' 0) w){i^ K ' 0, wi +m)} 

<Ce- cn z sup W?(x^' 0) ,£',P + /V 



(H-6) 

where Mi — number of points in x^ li ^ which inter- 
act with g-" 5 '°^(t) and p n( d = V n( (T)dT: recall that 
Vn s (t) = C jji ii^ 2 , see Eq.(6.10) (using the better esti- 
mate Eq.(6.11), V n < C vi wk , valid for all £, would lead 
to the same end result because (d > 1). Actually e~ cn ^ 
could be replaced by exp — (c2™«/ 2 ) as in Eq.(8.6)-(8.8). 

Consider first the case of p large, say p > p ni = 0(n^). 

Then W^x^; £, p) can be estimated by remarking that 
the argument leading to Eq.(6.9) remains unchanged if 

R(t) =P + Io V n ( ( T ) dT and R (t> S ) = R (t) + f s V n ( (r)dT 

are used instead of the corresponding R n( (t), R n( (t, s) 
(as long as p > 0). Then 



W{x { ^;^p + p ni ) <CW(x,p + 2p n( ) (11.7) 

as in the first of Eq.(6.8). 

Suppose po - p n6 > gc(£/7v), i.e. if p > Cn^, then 

W(x{ a) ;t,Po)) <C'W(x,p +p ni )) <C"(p + Pn( )) d < 
CpQ-. hence only the values of (n^ — l)^ < p a < Cn^ are 
still to be examined. 

In this case, however, the bound W(xf , p + p ni ) < 

CW(x,p + 2p Il? ) involves quantities p, p n( with ra- 
tios bounded above and below by a constant, hence 

W{x [ t ] ; £, p) is bounded by W(x; ^p + Cn^)< C'p d . 
Conclusion: there is C > 0, depending only on £q and 

for all p > gc(t/r v ), t < 6 it is Wix^-^p) < Cp d . 



D. Appendix: Quasi invariance 

A probability distribution p on a, piecewise regular 
manifold M is quasi invariant for a flow x — > Stx gen- 
erated by a differential equation x = v x if e- A « < 
p(S- t dx)/p(dx) < e A ^ and A(f) < oo. 

Suppose given O > 0, a piecewise smooth surface E C 
M with unit normal vector n x and a "stopping time" 
x — > #(x) < 9 defined on E consider all points x e E 
which are reached /or i/ie /irsi irnie in positive time t < 
#(x) from data j/ ^ E. Call _E, the set of such points, the 
fw&e wif/i frase E and ceiling 6(x). 

The probability distribution ^ is quasi invariant with 
respect to E and to the stopping time x — » if 
it is absolutely continuous with respect to the volume 
measure, its density r(x) is continuous and e~ A < 
p(S- t dx)/p(dx) < e x for some A > and for all < 
t < 6(x): this is referred to by saying that p is quasi 



invariant with respect to the stopping time 9(x) on E: 
symbolically p is (E, #(x))-A-quasi invariant. 

Then the following Sinai's lemma, [11-13], holds: 

Lemma: If p is (Ti,9(x))-X-quasi invariant the integral 
of any non negative function f over the tube with base E 
and ceiling 9{x) can be bounded by 



r r r e ( x ) 

/ f(y)Kdy) < e A / / r(x) f(S- T x) v x ■ n x ds x dr, 

J E JS JO 



> e" 



Jn Jo 



r(x) f(S- T x)v x ■ n x ds x dr (11. 



The lemma can been used to reduce dynamical esti- 
mates to equilibrium estimates. 

Proof: Let the trajectory of a point y which reaches 
E within the stopping time at x e E be parameter- 
ized by the time r and let ds x be the surface element 
on E. Then the set of points into which the paral- 
lelepiped A with base ds x and height dr becomes a paral- 
lelepiped S- t A with base S t ds x and the same height dr. 



Therefore the measure of p(StA) is 
hence the integral of any positive function f(y) over the 
set E can be bounded above and below by the inte- 
gral of J s f(Stx)p(x)ds x dr if p(x)ds x dr is the mea- 
sure of A: the latter is r(x)v x ■ n x ds x dr. Therefore 
p(x) =r(x) v x -n x . 



E. Appendix: Proof of Eq.(9.12) 

The factor e~ c Af2y arises because of the entropy 
bound (i.e. from the phase space contraction estimate 
within the stopping time). Therefore it is sufficient to 
find a bound to the integral in Eq.(9.10) without the fac- 
tor e^"'*). 

Consider first the case of S^. By Eq.(9.11) if y € <S| 
then |y n | = and there is (g, q) € y with q e dC^. 

Remark that y is the configuration reached starting 
from an initial data x £ Re within a time < Tm,V;u(x) < 
tA n (x): hence Eq.(6.11) applies. By Eq.(6.11) w(y) < 
\q\ < vi C y/n so that 

ff)(y) 

Pax'(dy) / dtw(y) 

(11.9) 

/j 
p(dx) 1 
Zc^x) 

where p(dx) is the ^io-distribution of configurations x 
outside Q and 

(11.10) 
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The estimate of the r.h.s. of Eq.(11.9), as remarked, 
is an "equilibrium estimate" . By superstability, [3] , the 
configurational energy U(q\x) > bk 2 — b'k^, so that Ji is 
bounded by: 



Be -ft(Mg-t'fc { ) IQI ' \ dC s\ r^L^k ( (n n) 
(k^ — 1)! fym 



while J jj,{dx) 



< 1 because ZcAx) > 1: and the 



bound can be summed over k^. Thus the contribution 
from Si to Eq.(11.9) is bounded by 



C' e GM 2 V /^ e -b[(logn) A 9A (?/r v )] 2 



(11.12) 



with C, b suitable positive constants. Since A > 1/2, this 
is summable over £ and yields the part of the Eq.(9.8) 
coming from the integration over <S>| . 

Let, next, y € <S| and let (g,g) as in (9.11). The func- 
tion w is 



it; 



|rf£(g,g)M| 
\gradE(q,q)\ 



(11.13) 



|<iE(g, g)/cit| < C Igln 1 / 2 because dE/dt is the work on 
the particle (g, g) done by the pair interactions (exclud- 
ing the wall forces). It is then bounded proportion- 
ally to the number of particles which can interact with 
(g, g), which, by theorem 4, is bounded proportionally to 
n 1 / 2 (as the total configuration is in I]'). On the other 
hand, |gradi?(g, g)| = y/\dip(q)\ 2 + m\q\ 2 > \/m\q\ hence 
w < Cn 1 / 2 again by Eq.(6.11) and the remark preceding 
Eq.(11.9). 

Then, analogously to (11.9), the integral under con- 
sideration is bounded by Ce c M v y/n (C, C are suitable 
constants functions of E) times an equilibrium integral 
J p(dx) Z c\x) w hh J 2 defined by: 



E 



■k — l w n$ — i 



dq 2 . . . dqkdq 2 . . .dqk 



e-h* area({ J B(g,g)= J B £ }) 



e -0jK(q 2 ,--,qk) 



(11.14) 



where the &iea({E(q, g) = E^}) is the area of the surface 
{(g, g) : E(q, q) = E ( } in R 2d (the is defined in (9.11)). 
Then J 2 is bounded by 



Ett-^ttH^, 



2tt s d\(*-i) 



fc v- • ^' 

(11.15) 

so that, suitably redefining C,C (functions of E), the 
contribution from <S 2 is bounded by 

C"e CM2v V^e-^ lo ^ X ^(i/r,)] 2 (1L16) 
and Eq.(9.8) follows from Eq.(11.12) and (11.16). 



F. Appendix: Regularized thermostatted dynamics 

Consider N particles in Uj ttj n A with a configuration 
of immobile particles outside A. The analysis in [7] can be 
followed and the solution of the equations of motion can 
be defined on the set T + consisting of the configurations 
x in which 1 of the particles is at £ on the boundary 
<9A, where clastic collisions take place, with normal speed 
q ■ n(£) > 0. The time evolution makes sense until the 
time t + (x) of next collision; it can then be continued 
after the elastic collision because, apart from a set of 
zero volume, the normal speed of the collision can be 
assumed ^ 0, until the time t*(x) > 0, if any, in which 
the total kinetic energy in one of the containers Clj, j > 0, 
vanishes. 

Remark that even in the cases in which the kinetic 
energy K m i n (x) d = min^o Kj^S^' 1 ^ x) vanishes as t — > 

t*{x) < +oo, the limit as t -» t*(x) of K jA (S^ A,1) x) and 
of all speeds and positions will exist (because the accel- 
erations otjCiji are bounded by max \dcp\N 2 uniformly in 
Kmin, using Schwartz' inequality). 

Hence a map T between T + into itself, mapping one 
collision to the outcome of the next, is defined for almost 
all points of T+ , [7] , unless the point x evolves into one 
with K m in(a;) = 0. 

Restricting attention to the points of x € T + whose ki- 
netic energies in any thermostat do not vanish for < t < 
t + (T p x), as in [7], for p < p* the dynamics S^'^x is well 
defined up to the time 6*{x) = Y?p t+(TP X ). The value 
of p* is p* = +oo unless for some value p* the particles 
of T p x grind to a halt before the next collision (which 
would, therefore, remain undefined since the equations of 
motion become signular). In the latter case a time t\(x) 
is defined signaling the moment in which the singularity 
occurs (an event not considered in the quoted reference 
because in the Hamiltonian equations considered there 
was no singularity of this kind). 

Until the time t\(x) > 9*(x) the dynamics exists and 
the only question relevant for us is whether 9*{x) < 
min(iA(a;), O): this would mean that there are infinitely 
many collisions with the walls and [0,9* (x)] would be- 
come the natural time of existence of the evolution rather 
than the smallest between tA{x) and 9, as used in this 
paper. 

Suppose that 9*(x) < 9 and that K j>A (S^ A,1) x) > \ 
for t < 9*(x): call such points x. Then the volume 
contraction of the distribution JIq obtained by condition- 
ing /Jo to the particles outside A is bounded uniformly in 
all subintervals of [0,9* (x)): hence a set A in gen- 
erates a "tube" A* = U xE a ^o<t<t + (x) S^'^x and the 
volume ~fi (T n A*) > A7J (A*) where A is a lower bound 
on the 7Z -volume contraction over any time interval in 
[0, 9] over which the dynamics is defined. The quantity 
A is bounded for all n such that YT P =o t + (TP X ) < 9*{x). 

Therefore the sets T k A* cannot be disjoint for all k 
unless ^2 k T + {T k x) > t\(x): i.e. this remark takes the 
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place of Poincare's recurrence theorem used in [7] and 
allows us to conclude that until the thermostats kinetic 
energies are all positive the regularized dynamics exists 
and the elastic collisions with the boundary of A cannot 
accumulate in time. 



This means that the evolution proceeds until the first 
time t^(x) (if any) when some of the Kaj vanishes. 
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